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ABSTRACT 


Preliminary Turbine Design for Thermal Cycle Feasibility 
David A. Potter 

Submitted to the Department of Ucean Engineering and the 
Department of Mechanical Engineering on 10 May 1974 in partial 
fulfillment of the requirements for the degree of Master of 
Science in Naval Architecture and Marine Engineering and 
Master of Science in Mechanical Engineering. 

Mathematical models of an axial flow turbine and a radial 
inflow turbine were developed. Utilizing these models, a com- 
puter program was written such that given a thermal cycle and 
@ minimum of input variables, the necessary turbine design cal- 
culations are performed. An extensive output is provided such 
thut the designer may utilize this output in order to evaluate 
the practical feasibility of the given thermal cycle. 

The Zener sea cycle is used as a test case with axial and 
radial designs examined for feasibility. 

The computer program is written in the Fortran IV language 
specifically to be run on the Interdata 70 computer. 
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NOTATION 


Blade cross sectional hub area 
Flow area 

Annulus 

Blade chord length 


opecific heat at constant 
pressure 


Specific heat at constant 
volume 


Hydraulic diameter 
Gravitational constant 
Enthalpy 


Moment of Inertia in 
tangential direction 


Moment of Inertia in 
Axial direction 


Conversion factor 
tindece loss 
Stator loss 

kotor loss 
Clearance loss 
Mass flow rate 
Mach Number 


Bending moment in tangential 
direction 


Bending moment in axial 
direction 


UNITS 
ft? 
ft" 
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rt- 


BTU/°F 
BTU/°F 


ee 
B22l74 fie rem, wpt were 
BTU/1bm 


rt 


et 


778 ft-lbf/BTU 
BTU/1bm 
BTU/1bm 
BTU/1bm 
BTU/1bm 
lbm/sec 
Non-Dim 


ft-lbf 


ft-lbf 
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Number of axial stages 
Specific speed 

Throat opening 
Pressure 

Volumetric flow rate 
Universal gas constant 
Reaction 

Radius 

Reynold's Number 

Shaft rotative speed 
Spacing 

Entropy 

Temperature 

Blade height 

Rotor tangential speed 
True fluid speed 
Relative fluid speed 
Axial distance 


True fluid flow angle from 
axial 


Relative fluid flow angle 
from axial 


Ratiomot Cc /C 
p Vv 


Primary loss factor 


Primary loss factor corrected 


for aspect ratio 


Non-Dim 
Non-Dim 
feet 
lbf/in® 
(a ae 


CH 


BYU/1 bra 
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feet 
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Rev/min 
feet 


BTU/1bm © 


R 
Deg Rankine 
feet 
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ft/sec 
ft/sec 
feet 


Deg 
Deg 
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Non-Dim 


Non-Dim 
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aS ee ee factor corrected Non-Dim 
or AR & Re 

3 Slip factor Non—Dim 
y) Efficiency Non-Dim 
© Tangential angle Deg 

a Stage loading coefficient Non-Dim 
pel Static viscosity lbm/hr ft 
\ Kinematic viscosity Btw ee 
A Differential function Non-Dim 
oC solidity Non-Dim 
a. Centrifugal stress igen 
YT Constant 3.1415927 
p Density lbm/ ft 
d Flow coefficient Non-Dim 
YY Work coetficient Non—Dim 
ee) Angular rotative speed Radians/sec 
SUBSCRIPTS 

O Stagnation conditions 

ol Machine inlet stagnation condition 

02 Machine outlet stagnation condition 

o2s Isentropic outlet stagnation condition 

1 Individual flow path 

t Total flow paths 

v) Tangential direction 


X Axial or meridional direction 
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mat. 
max 
Axial 
1 
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Radial 


& WW NW 


Mean 

Hub 

Tip 

Fluid 

Total to total 
Total to static 
Material 


Maximum 


Rotor inlet 


Rotor outlet 


Stator inlet 
Stator outlet 
Rotor inlet 


Rotor outlet 
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INTRODUCTION 


The thermal cycle is the basic tool in use today to con- 
vert the energy contained in a source or fuel into useful 
power. Cycles using geothermal energy, cycles powered by 
solar energy, and cycles powered by the thermocline existing 
in the oceans, are just a few of the proposed cycles being 
studied as possible sources of power for the future. asi the 
search for new and plentiful energy sources in the years to 
come, new and different thermal cycles will be proposed. 

The evaluation of a given thermal cycle as to its prac- 
tical feasibility involves a detailed study in which not only 
the energy source, but the hardware needed to convert that 
energy into useful power, are investigated. In such a study, 
the engineer must concern himself with three main areas: con- 
verting the energy source to heat, transferring that heat 
energy to a working medium, and then utilizing the increased 
energy of the working medium to produce useful work or power. 
It is with this last facet of the problem that this thesis 1s 


concerned. 
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BACKGROUND 
1.1 Thermal Cycles 

Thermal cycles come in all sizes and shapes. They are 
used not only to convert heat energy to useful work, but to 
transfer heat energy from one location to another, as ina 
refrigeration unit. The difference between the two cycles, 
other than hardware, lies in the fact that in a work-orienta- 
ted thermal cycle energy is supplied in the form of heat or 
fuel and work is the end result; while in a heat transfer 
or refrigeration cycle, work is supplied as an input and the 
end result is the transfer of heat energy. While both types 
of cycles operate under the same principles, it is the former 
type which is of great concern in the search for means to 
supply society's growing energy demands, and which will be 
the only one considered from hereon. 

Thermal cycles of the heat to work type can be catego- 
rized as either open or closed cycles. In the open cycle, 
the energy of the source (fuel) is transferred to a working 
medium in the form of heat, whereupon this increased energy 
of the working medium is utilized to produce useful work. 

At the conclusion of the open cycle, the working medium is 
exhausted from the machine. In a closed cycle the process 
1s Wee satially the same; however, the working medium is re- 
tained in the cycle after producing work and is reintroduced 
to the source where it again increases its energy state with 


the addition of heat. There are, of course, many variations 
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and combinations of the above cycles using many different 
working mediums at many WSU suet temperatures and pressures. 
1.2 Turbomachinery 

One problem common to all thermal cycles, open or closed, 
is that of finding an efficient means to convert the energy 
of the working medium into useful work. Without considering 
unproven schemes, such as magnetohydrodynamics, there are 
two basic classes of machines in use today with which this 
can be accomplished: machines which extract the energy from 
the working medium as a change in its potential energy, and 
those machines which extract the energy as a change in the 
kinetic energy of the working medium. Without discussing 
the relative merits of either type of machine, it is suffi- 
cient to state that the second class of machine, the turbo-—- 
machine, is a viable alternative and one which should be 
considered in any thermal cycle feasibility study. 

Turbomachines, or more specifically, turbines used in 
thermal cycles today, can be divided into two categories: 
those in which the working fluid flows through the machine 
in a nearly constant axial direction parallel to the center- 
line of the machine; and those in which the working fluid 
flows in a direction radial to the centerline of the machine, 
either inward or outward. ‘The choice as to which type of 
machine is best for a certain set of conditions is a diffi- 


cult one, depending on a large number of variables external 
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to the cycle itself, such as weight, cost, and operating 
load. This will be discussed more fully in following 
sections. 
1.3 The Design Sequence 
In determining the technical feasibility and practical 
value of a given cycle, at the preliminary design stage, the 
engineer must follow an iterative process. One logical 
sequence would be to: 
1) Assume a-working medium and a set of states for 
the medium (most commonly the energy extraction 
process is done in a gaseous state). 
2) Assume certain hardware restraints and a flow 
rate for the working fluid. 
3) Design the required hardware (if not available 
commercially). 
4) Evaluate the necessary hardware and if not 
Satisfactory return to step l. 
This process is repeated several times until the hardware 
conforms to the engineer's definition of feasibility, or until 
some prechosen parameters are optimized. The problem here 
lies in the fact that the calculations necessary to design the 
hardware, although relatively straightforward, are numerous 
and time-consuming. For a given set of cycle conditions, 
there are often many options open to the designer, all of which 
will result in different hardware and all of which will meet 


the design criteria. The choices between axial or radial 
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turbines, and between turbines or reciprocating machines, 
are good examples of this. There is a necd, therefore, for 
the designer or engineer to quickly translate his cycle re- 
quirements into specific hardware sizes and characteristics 
in order to logically chvuose between the options available. 
1.4 Optimization 

Feasibility and optimization are two words whose’ 
definitions can be different for every different cycle con- 
sidered. Ina given situation it may be necessary to max— 
imize efficiency, while imu another situation, the design 
may be constrained such that weight or volume or cost must 
be minimized. As for turbines, in some applications an 
axial flow unit might be chosen for its high efficiency, 
while in others a radial machine with its lower fabrication 
cost might be best. Feasibility also has different meanings 
in most different contexts. In one instance it might mean 
technically possible, while in another it might mean finan- 
cially profitable. In order to auswer these questions, the 
engineer or designer must first define feasibility and then 
decide upon which, if any, parameters are to be optimized 
within the feasibility constraints. Having done this, he 
may then cvaluate the hardware designs and make a logical 


choice betwcen competing designs. 


1.5 Purpose 


Within this framework then, the purpose of this thesis 


is to bring together the various concepts of turbine dcsign; 
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to develop a mathematical model; and to computerize that 
model in such a way that, given a minimal amount of in- 
formation about the cycle and the machine constraints, a 
set of turbine hardware parameters will be provided that 
meet the requirements of the preliminary design stage. 
These may then be used in the design sequence mentioned 
previously to determine feasibility. No attempt willbe 
made to optimize within the computer program, but rather 
to provide information from which a suitable optimization 
criteria can be evaluated. for this reason no attempt will 
be made to choose between an axial or a radial turbine. 
Both will be developed simultaneously and information will 
be provided on both types. The program will be capable of 
handling multiple inputs for the same cycle in one computer 
Fun < 

It is hoped that by utilizing this program the design 
sequence may be shortened considerably. This will be true 
especially if use is made of other computer programs to 
calculate heat exchanger parameters. In doing so, the 
designer will be able to study more completely the possibi- 
lities of a given cycle at the preliminary design stage, 


and thus develop a more suitable final design. 
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DEVELOPMENT OF THE MATHEMATICAL MODEL 
2.1 General 

A design philosophy 18S a necessary ingredient to any 
worthwhile design. This philosophy may not always be ex- 
plicitly stated; however, it is always present in a good 
design. It is a means whereby the designer is able to make 
logical and, of most importance, consistent decisions during 
the course of the design sequence. This is especially true 
at the preliminary design stage where decisions are made 
which constrain the course of the design throughout. 

The philosophy followed in compounding the mathematical 
models described in the following sections and in the actual 
programming itself can be stated as accurate and simple. 
Both of these qualities apply to any work done at a prelimi- 
nary design stage. Accuracy 1S required because the designs 
have now passed beyond the back of the envelope stage; and, 
although detailed design is not required yet,.the results 
must be accurate enough so that detailed design can begin by 
using them as a basis. Simplicity 1s a must because at the 
preliminary stage, neither the time nor the money are usually 
available to carry out the extensive calculations of a de- 
tailed design. That is not the purpose of the preliminary 
design. One last idea is also followed in the handling of 
certain relationships; that in the absence of accurate and 
Simple relationships to cover a certain situation, traditional 


and reasonable approximations are assumed. With the lack of 
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specific data in certain areas this borders on au educated 
guess, but then again this is what the designer is paid for. 

As to the development of the models themselves, they 
can be divided into three main areas: representation of the 
working medium itself, representation of the axial turbine, 
and representation of the radial turbine. They will be ex- 
plained iu that order. 

It is assumed that the machines described by the models 
can be built of conventional materials and that the axial 
blading has no twist and is set ina purely radial direction. 
It is further assumed that internally cooled blades or rotors 
are not utilized. 

All computations are in English units with the notation 
as explained on pages 8 to ll. Calculations are carried 
through to a precision consistent with the accuracy of the 


Interdata 70 computer (six significant figures). 


2.2 The Working Medium 
2.2.1 Equation of State 


Two assumptions concerning the working medium are basic 
but should be stated for clarity. It is assumed that the 
workiug medium is a fluid and that it enters the machine in 
a gaseous state; this precludes the design of hydraulic tur- 
bines. The second assumption is that the inlet and exit 
states of the fluid are known and that the designer has the 
Breer ti cls of the fluid at these states available as input. 


In order to accurately represent the fluid within the 
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machine, an equation of state for the fluid had to be picked. 
The perfect gas model was chosen over the other equations of 
state for its flexibility, its accuracy, its wide acceptance, 
and its compatibility with the standard one dimensional com- 
pressible flow relationships. It is therefore assumed that 
in the region under consideration, the fluid can be repre- 
sented by the perfect gas model and that the effects of 
velocity can be represented by the standard one dimensional 
compressible flow relationships. 

To make this assumption more realistic, it is necessary 
to calculate the different constants in the region under con- 
sideration rather than use traditional or tabulated values. 
One method to achieve this is to use data for an isentropic 
expansion in the region and to calculate the constants from 
this. Since the designer already has calculated the point 
in order to determine his outlet conditions, it was decided 
to use an isentropic expansion across the machine for the 


reference points. Thus from points 1 to 2s, 


5 = (ooo) /(To9577T 91) 


Using the Gibbs equation and the perfect gas law, the 
following commonly used relationship for an isentropic 


expansion may be derived: 
ml 3*/ (x-1) 
(Poa s/ Poi) - ce ey 


where WESC IO 
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Solving this equation for *%: 


is In (Po 5 /Po))/(inl(P5./P)) 7 in(T)5./T),)) 
Thus C= om) a | 


and He= “C= CC. 
P V 


2.2.2 Flow Kate 

In describing a particular turbomachine, the mass flow 
rate of the machine has a large effect on its design. For a 
particular cycle, the mass flow rate to achieve a given 
output can be easily calculated using the steady flow energy 
equation. It may not be advantageous, however, to achieve 
this output on only one machine, but rather to utilize several 
Iachines, connected in parallel, to achieve the given output. 
The mass flow rate through a given machine, with no allow- 
ance for leakage, wili be: 


Mo gayi dual tot / # units in parallel 


2.2.3 Compressible Flow Relations 

The two parameters necessitating the use of the one 
dimensional flow relations are the calculation of Mach number 
and the calculation of flow area. Using these relationships, 
and given the stagnation or total temperature and the 
velocity of a fluid at a point, the Mach nuinber and flow 
area may be easily calculated. From the definition of total 


temperature and Mach number: 


2 | 
pot sy /2¢ 56, 





Aap 


as 
i ee (sg kT) * 


in His common form: 


iis (T,)?/PoAg = M (ag, /R)?/ (1+ (8-1/2 )M2) P42) /2 2) 


Thus A 


1 
f (mh (T)?°/P.) x £( 3g ,R,M) 


‘AL. 
Bs 


i 


1 
o 2 e 
(nm (T) La) x (flow function) 


where: 


Flow function = (1+ (x1 /2)M2) PHD) /2 (098-1) jy (ag /it)? 


Bee Axial Turbine 
2.3.1 General 

The first step in the design sequence involves deciding 
the general configuration of the machine. Is the hub radius, 
the mean radius, the blade tip radius or some combination of 
these to be maintained constant throughout the machine? 
Wilson (1) recommends that the mean radius be maintained 
constant, while Carmichael (2) recommends hub radius be kept 
Monstant at the preliminary design stage. Constant mean 
radius was chosen for the reason that the effect on turbine 
performance between the two is small, and calculations at a 
luter stage will be simplified with a constant mean radius. 
Following this assumption, mean conditions are calculated 
first. As the velocity triangles and other blade parameters 
vary by stage throughout the machine and as the extremes are 
represented by the first and last stages, only the first and 


last stages will be calculated in detail with the other 
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stages assumed to lie in between these values. Finally, 
machine parameters which depend upon the previous calcu- 
lation, such as total efficiency and machine length, will 
be made. It is also assumed at this point that the machine 
is a full admission turbine operated under a constant load; 
both assumptions are consistent with stationary power plant 
operation. 

2.3.2 Reguired Vutput 

By a process of elimination, the list of quantities 
which is sufficient to describe the turbine at the prelimi- 
nary design stage was narrowed to the following. Many other 
quantities could have been included; however, all major para- 
meters are included, and those omitted may simply be calcu- 


lated from those given. 


1) First and Last stage Tip Diameters 

2) " Blade Height 

3) : Aspect Radios 

4) _ Number of Blades 

5) Hy Efficiency ( Total to total) 
6) iL Velocity Triangles 

7) Machine Critical Mach number 

8) u Efficiency ( Total to total) 
9) : Critical Stresses 


10) Length 
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In describing the velocity triangles for a given stage, 
three locations for both the stator and rotor entrance tri- 
angles are described: hub, mean and tip. 

a. se) input 

in listing the factors which are generally used to des- 
cribe a turbine stage, it was found that they could be cate- 
gorized in three groups: those which describe the mass flow 
rate; those which describe the work accomplished per pound 
of fluid per stage; and those which describe the performance 
of a given stage. Each facet is necessary for the preliminary 
design. 

The mass flow rate for the machine, and thus each stage, 
(since the stages are assumed to be in series with no bleed- 
offs) can be described by the total mass flow rate for the 
cycle and the number of turbines in parallel, as described 
mm ssection 2.2.2. 

The work per pound of fluid per stage can be simply re- 
presented also. Making the assumption that the work accom- 
plished per stage is equal for all stages (2), the individual 
stage work can be represented by the total work or enthalpy 
drop for the machine and the number of stages. 

The parameters used to describe stage performance are 
not so straightforward. They can be classed as either pri- 
mary or secondary, depending upon the strength of their in- 
eters upon performance. Primary factors are generally 


associated with the geometry of the velocity triangles, and 
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if the deviation is assumed negligible, the geometry of the 
blade itself. Secondary factors have a lesser effect on 
efficiency and are related to the likes of clearance, aspect 
ratio, and other mainly mechanical factors. In order to main- 
tain a simple machine only primary factors will be specified 
as an input, with a logic developed for the determination of 
reasonable secondary factors based on the given primary 
factors. 

As primary factors, Carmichael (2) recommends reaction, 


, and inlet flow angle, x Craig and Cox (3) prefer to 


ya 
use the stage loading parameter, Y, also called the work 
coefficient, and the flow coefficient, ®. A is another widely 


used parameter (4). Comparing the definitions of these 


parameters, 


R= 1 - (Vg) + Vg9)/2U 
Tan (%,) = V/V, 

y 
® 


A. 


va 
AUV,/U 


II 


vee 


U/ (Voy 7 Vga) 


it can be seen that there are only four independent variables. 
For reasons explained in a later section, one .of these can be 
eliminated; thus only three independent variables will des- 


cribe the shape of the velocity triangles. 
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The three variables chosen for use here are R,U ; and % - 
Reaction was chosen for its compatibility with the Euler work 
equation which will be used later to solve for a mean radius. 
Given this mean radius and an KPM, the second variable, U, . 


can be specified. The last variable, & was chosen as it 


1? 
would heave to be calculated anyway during the course of the 
design. 
The input for the axial turbine can be summarized as: 
1) Mass flow rate 
2) Number of units in parallel 
3) Number of stages per unit 
4) Shaft revolutions per minute 
5) Reaction at mean conditions 


6) x 


l at mean conditions 
7) Cycle conditions of section 2.2.1. 

2.3.4 Mean Conditions 

The starting point for the mean conditions is with the 
Simultaneous solution of the Euler work equation and the 
reaction equation. It is assumed that the fluid enters the 
stator only in an axial direction, with no tangential compo- 
nent. This is a common assumption (4) and is made for several 
reasons. One argument is that the turbine works by extracting 
energy from the fluid; thus any energy left in the fluid at 


rotor exit (stator inlet) in the form of vorticity (tangential 


velocity) is unavailable for extraction, and thus contributes 
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to a loss. Using this assumption, the Euler equation applied 
to mean conditions across the rotor becomes 


2 Le 7 © © Voi m/ 8% 


n 


where Aho = hoy 862 


The reaction equation, also applied at mean conditions across 


the rotor 1s simplified to 


Che Vo) ,/24T,, + 


The simultaneous solution of these equations for I and Vom 
with W=NRPM/30 gives 


ap 
r= (ah g,J450/r RPM n (1-Te) ) ? 


Wein = rYRPM(1- R)/15. 


Substituting in the appropriate constants, 


1 
as ei 2 /p py 
a = 1068.247(ah /n(1 R)) /RPM 


and 


1 
= - o3 
Me ne 223.7331 (4h (1 n)/n) : 


In order to prevent the solution from expanding at Ro=1, 
must be limited to less than 0.99. 
The rotor speed can be calculated as 


U, = U, = U = YT RPM r / 30. 


Using the trignometric relationships, Figure 1, and the 
eo, which was given as an input, the rest of the velocity 


triangle for the rotor inlet can be calculated. 
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Axial Turbine Velocity Triangles 


Rotor Inlet 





Stator Inlet 





Figure 1 
Yom = Yorn/Sin(,,) 
Vi = Voin/ Tantes,) 
A, = ArcTan(Tan(%_) (1-1/(2-2R.) )) 
= V /Cos(A, ). 


Making the assumption that the axial velocity is constant 
throughout the machine as was stated before, the tangential 
velocity at rotor outlet is zero; thus the velocity triangle 
for the stator inlet (rotor outlet) can be calculated: 


P50 = ArcTan(Tan(& _)/2(1-R )) 





oo 


Wor .: Voim/ Vos (A,,,) Tan (04 ,) 
Om 7 oon — 

Yu a ee - Vs 

on = Vio 


The total gas path deflection for the rotor is the sum 


of the inlet and outlet angles CA +f,_)3 while for the 


of 3 
lm 


2.3.5 Stage Conditions 


stator it is just 


‘as the inlet angle is zero. 


In solving for stage conditions it can be assumed that 


the extreme values are represented by the first and last 


stages, with all other stages falling in between. 


Variance is assumed linear, the 
can be calculated for any stage 
This assumption also allows the 


efficiency as the mean of first 


Based on these assumptions, the 


If the 
dimensions and velocities 
with reasonable accuracy. 
calculation of machine 

and last stage efficiencies. 


following calculations can 


be made for the first and last stages. 


Using the one dimensional compressible flow relationships 


described in section 2.2.3, the 
@a given velocity and stagnation 
and annulus area are related as 


Area = Bales 
ann 


flo 


flow area can be calculated for 
temperature. The flow area 


shown in Figure 2. 


/Cos(&). 


Usiug the definition of annulus area and the mean radius, 


the blade height may be found. 
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Axial Turbine Flow Area 


Machine 


Axis . 
a> 





Figure 2 


Z 2 
Area =\f(r, ae ) 


=Yi(r <r, ir +r, ) 
=‘Nft2r 
m 
Blade height = t = Area Pa 39. Say 
ann m 


where r= (xr tr, )/2 


The tip and hub radii are found by respectively adding or 
subtracting one half of the blade height to the mean radius. 
The next step involves making the assumption that a 


streamline of fluid has no appreciable curvature in the radial 
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direction as it moves through the machine. This is not com- 
pletely correct, but for small wall curvatures it can be con- 
sidered accurate (1). This assumption means that a segment 


of fluid, as shown in Figure 3, must be in radial equilibrium. 
RADIAL EQULLIBRLUM FORCES 


p + dp dr 
dr 


CW 


/ 


f yi 

/ / as dr 
o i 

i ye f 

Pp 

deg 


V, 


Xxis of machine into page 


Figure 3 
Summuation of the forces acting on this segment gives the 
equation for this condition of radial equilibrium. 


2 
ha lity = ge, 


/¢ dr 
Using the Gibbs equation and the definition of total 


enthalpy for an adiabatic condition, the radial equilibrium 
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equation can be rewritten as 


dh -—- Tds = 1 dV ‘ + ] d(r°V ae 
O x 2 3) 


dr dr 28 dr 2g or dr 











Making the assumption that the enthalpy drop is uniformly 
distributed over the blade radially, and that entropy is not a 


function of radius, the result is 


oY 
dV é = =) Are iy) 
= A 7 | ae 
dr r ar 


At this point one of several distributions for Nee or Vo 
may be chosen; however, having made the assumption previously 
that axial velocity 1s constant throughout the machine, the 
above reduces to 

i Vo = constant. 

This constant may be calculated from the mean conditions 
and, using the above relationship, the tangential velocities 
at hub and tip can be calculated. The blade speed is easily 
calculated at hub and tip as shown previously; and, using 
trignometric relationships, the complete velocity triangles 
for rotor and stator are easily calculated in the manner 
illustrated in section 2.3.3. 

The choice of aspect ratio (blade height/blade chord) 
must now be made. As this is a secondary variable, section 
2.3.3, a means must be derived for a logical decision. A 
study of duta presented by Forrester (5) shows that the 


secondary losses attributed to aspect ratio effects decrease 


with an increasing aspect ratio; thus a large aspect ratio 





be, 
would be advantageous from an Rrerercnee standpoint. On the 
other hand, Horlock (6) states that the bending stresses on 
the blade may place a limit on the blade chord length and 
thus limit the feasible aspect ratio. Any aspect ratio 
chosen should therefore be as large as possible without ex- 
ceeding a predetermined stress limit. Carmichael recommends 
a limit of one ton per square inch which was used in the 
following derivation. 

In order to simplify the calculations necessary, the 
fluid was modeled as incompressible and the blade shape was 
Simplified so that the resulting criteria could be applied to 
the general case. The maximum stress was assumed to occur in 


the rotor at the root of the blade and located on the 


Blade Root Stresses 





Figure 4 
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upstream side. By considering the attrtce in tangential momen- 
tum across the blade, the bending moment in the tangential 
direction (Figure 4) can be, shown to be 


z 
— t = r 
Mg =f “AeVy2 tr Meg) Cbs 
2 ees 
h Z 
where z = number of blades per stage. 
Integrating this equation, and substituting the Luler 
work equation and the definition of blade speed gives 
: 2 
iv = = 
My = floh J(r, -r,)"60 S 
ee 


i . 


-r,*) (21r)° RPM r, n 


h 
Simplifying and substituting in the constants yields 


M.~ = 591.209 hah | 


Ph 
9 ) S,/RPM rT, 2 


*t *h h 


Since Tan (2) = Forces / Forces 


6 direction x direction 


the bending moment in the axial direction can be described by 
M. = My Tan (4). 


For a particular design let Kl be a constant representing 


Kl = ah | is a5) 591.209/RPM a a 
then Mo = Kl S35 
M, = Kl Sy Tan (A). 


For the blade shape shown in Figure 4, the moments of 
inertia about the respective neutral axes can be shown to be 


Tee s,°/69.429 


@! 


3 
eS ae 
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For stress at the point A, 


2 
@/ 1s 


30.86/c Sj) 


2 
x/Tz 36/S)c¢ e 
The bending stress is then 


o-, = Ki(Tan(4)36/c~ + 30.86/c S,). 
Using Zweifel'‘s criteria to relate S and c (4 


0.8 = 2 Cos(4,)(Tan(Z, ) + Tan(4,, )) s/c 


substituting in the bending stress equation and solving for 
chord, using the one ton per square inch stress limit men- 


tioned previously, 


nl— 


c = (K2mif) 
where K2 = Kl 32.174/(2000) (144) 


fe 36 Tan(Q, ,,) a (eee Cos” (G5, ) X 
(Tan(4), ) + Tan(4,,)). 


This approximation for chord is based on a stress limit 
caused by gas bending forces; however, to account for the 
fact that other design conditions may impose a lower limit 
than this criteria, a maximum limit to aspect ratio of 10 
will be imposed. 

Zweifel's criteria can now be used to determine spacing 
for the rotor from this calculated chord. AsSuming the same 
chord for the stator, spacing can be calculated there also. 
As shown in Figure 5, the throat opening can be figured for 


Meo the rotor and stator. 
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Throat Opening 


Axis of machine 


W 
—\ O Nee 
Ss 


ay 


Figure 5 
0, = 5) Cos(_) 
0, = 8, Cos(f,_) 


Knowing the spacing, the number of blades for rotor or 
stator is simply 
i arr, / Se 
yeine the hydraulic dianeter of the throat opening, 
again for both rotor and stator, 
jb) ER OR ae) 


the Reynolds number for both the rotor and stator can be 
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calculated as 


ke 


2 = Dy Wo/% 


lie, 


D, V44/%- 
Bao .o0 bificiency 

The determination of efficiency for axial turbines is a 
much-debated subject. Craig & Cox 0 Sina tien) , and Ainley 
and Matheson (5), all correlate their data using the work 
and flow coefficients mentioned previously; while Soderberg 
(4) correlates using gas deflection angle. In either case 
the primary loss is related to blade profile loss. Secondary 


losses are attributed to many factors, some of them being: 


1) Clearance loss 

2) heynolds number effects 
3) Aspect ratio 

4) Guide loss 

5) Runner loss 

6) Guide yvland loss 

7) Lacing wire loss 

8) Balance hole loss 

9) Wetness loss 

10) Dise windage loss 
11) Partial admission loss 
2 ean lias a lle los 


igjmeCava ty sloss. 

While including all of these factors would give a most 
accurate efficiency prediction, the accuracy required in a 
preliminary design doesn't require the inclusion of them all. 

Soderberg's correlation was chosen in that it provides 
accuracy consistent with the intent of the program; yet it 
1S simple enough to be easily programmed on a computer. This 


correlation takes into account profile loss, cleurance loss, 
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aspect ratio, and a Reynolds number effect. These losses are 
equated to a fraction of the kinetic energy available to the 
turbine and a total to total etficiency calculated. 

The method consists of finding a primary loss factor, 
based on gas deflection angie for both stator and rotor. 
This loss factor is modified for aspect ratio and Keynolds 
number effects, and then an efficiency calculated. The 
efficiency is then modified by a factor which takes into 


account the clearance losses. 


Po 
yy 


f (Gas Deflection) 


ea deg cafe oy tt) ed 
5 


(10?) 
(Oa ta (Re ) 





The above are calculated separately for both stator and rotor. 
The efficiency is then: 


ny = (1+ /,Stator V5“ {Rotor Woo) 


eu Vea! 


N.4 (Stage) = (Area Annulus - Tip Clearance Area). 
5 Area Annulus 


The prinary loss coefficient is one which can be pro- 
grammed easily in a curve-fitting subprogram. The tip clear- 
ance area was Calculated based on a tip clearance of one per- 
cent of the tip radius. This gave a constant factor of 0.99. 
Assuming a linear distribution of efficiencies between the 
first and last stages means that the machine efficiency is 


the average efficiency. 
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Were riat stare. |. Laet.Stere)/2- 


2.3.7 Miscellaneous Machine Parameters 

In order to more easily picture the machine size, all 
pertinent radii are converted to diameters. 

A machine unit length is calculated using the average 
stage chord and the number of stages. 

The maximum blade centrifugal stress occurs in the last 
stage at the rotor blade root. This can be written as 


follows: 


If the blade is assumed to have no taper, AL ls constant. 
4 

Substituting in the constants and assuming the density of 

steel 1s approximately 0.29 Lbm/in? gives 


5-) = ©.001086 RPM“ c™ +. 
Cc mM 


2.3.3 Turbine Limitations 

Certain limitations must be placed on the mathematical 
model previously described, so that it won't produce turbine 
parameters for a machine that 1s either rye cally impossible 
or not consistent with the assumptions of the model. These 
limitations will be in the form of numerical checks in the 
program and will cause the calculations to be -terminated if 
they ure not satisfied. 


The first check is for excessive vas deflection. While 
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this in itself is not an impossible situation, only leading to 
a lower efficiency, the data available for the Soderberg 
correlation only accommodates a gas deflection of up to 40-3 
Since the efficiency decreases rapidly at angles greater than 
this, no harm is done in limiting gas deflection to 140°. 

Keaction is another parameter which is limited. The 
upper bound has been discussed earlier in section 2.3.3. A 
check must be made, however, to insure that at no point does 
reaction become negative. Physically, a negative reaction 
would mean that within a given stage, at some point, a pressure 
rise was occuring to be compensated for by a larger pressure 
drop at another point. The overall condition, of course, 
would be a pressure drop. In order to avoid this condition 
which brings about higher losses, a check should be made to 
insure non-negative reaction. Negative reaction can only 


occur with large V, and since it is assumed that r Vg = 


6 
constant, the check should be made at the smallest radius, 
the last stage hub. 

In solving tor a mean radius, in section 2.3.4, no con- 
tinuity check was made. without a check for continuity it 
would be possible, using the model described, to compute a 
Negative machine diameter. In order to avoid this and to 
allow for bearings, shafting, and blade attachment, a check 


Will be made to insure that the hub radius calculated is 


Breater than 50 percent of the tip radius. 





41 

The one dimensional compressible flow equations used 
throughout also pose another Ae re tot The fluid flow 
throughout the machine is assumed to be subsonic, and 
the design has proceeded such that no consideration has been 
given to the shock effects caused by supersonic flow. To 
insure subsonic flow, a check will be made on the two largest 
velocities to insure a Mach erin less than one. The 
velocities are the last stage hub velocity at rotor inlet and 
the last stage tip relative velocity at rotor outlet. 

In addition to these above limitations, there are physical 
limitations imposed by the computer programming itself. These 
are listed in the appropriate appendices. 

2.4 lkadial Turbine 
2.4.1 General 

The term radial turbine is used to describe a wide range 
of turbine configurations. In the general sense, it describes 
a Meio niachine where the fluid enters the rotor with no appre- 
cliable axial component of velocity. Whether the fluid flows 
radially inward or radially outward and whether it leaves the 
rotor in a radial or axial direction depends upon the context 
of the discussion. The choice made here to limit the radial 
design to that of only one configuration was based on practi- 
cality. The 90 degree inward flow radial turbine, in which 
the fluid enters the rotor in a radial but leaves in an axial 


feeeetion, 1S the most common type of radial turbine (7). It 
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has the advantage of the strongest structural strength, also 
a factor which can be important at high rotational speeds. 
Only single stage turbines will be considered as these are 
the most common also. 

The methodology of designing the radial turbine can be 
summarized in the following four steps: rotor design, nozzle 
design, scroll design and efficiency correlation. The 
station designations are as shown in Figure 6. 

2.4.2 Output Kequired 

The parameters necessary to describe the turbine of 
Figure 6 can be narrowed down to the following list for 
preliminary design: 

Nozzle Inlet and Outlet diameters 
Nozzle blade height 

Rotor Inlet diameter 

Rotor Outlet diameters at hub and tip 
Maximum Scroll Area 

Unit length 

Specific speed 

Lfficiency 

Velocity triangles at 


A) Nozzle Inlet 
Bjmeotomeintet andweOutlet. 


OMONAULBRWN-E 
See” Sone” See” “eee” “Seem” See” See” “nee” “eee” 


The velocity triangles for rotor outlet will be described 
at hub, mean, and tip locations. 
oe4.3 Input 

The parameters necessary to describe a radial tiachine can 
be categorized in the same manner as the axial parameters 
(section 2.3.3): those which describe the mass flow rate, those 
which describe the work felcon pl iahode ee stuge, and those which 


describe the performance of a stage. 
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Radial Turbine Station Designations 


i 
stator —? jad 
3 





je Machine Axis 


Figure 6 


As it has been previously assumed, that the radial turbine 
design is for a single stave machine, the work accomplished per 
stage is the same as the total work to be accomplished by the 
evele perameters, section 2.2.3. 

The mass flow rate for this machine can also be described 
in a manner sililar to that for the axial machine and which 
is presented in section 2.2.2. ‘The variables necessary for 
this representation are the inass flow rate and the number of 
een nes operating in parallel along the fluid flow path. 

Performance estimation 1s more standardized for radial 
turbines than axial turbines. Through the use of a non-dimen- 


Sional specific speed parameter 1t is possible to correlate 
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efficiency with rotational speed, volumetric flow rate, and 
work accomplished, for turbines designed with the same basic 
assumptions and geometric relationships. (8) There are several 
specific speed parameters in use which will be discussed ina 
later section; however, for the immediate purpose, rotative 
speed is a necessary input. One last input remains and that 
relates to an assumption to be made in the succeeding’ section. 
Without explaining this assumption in detail, it is sufficient 
to say that the number of radial blades of the rotor should 
be included as an input. 

The input for the radial design can be summarized as 
follows: 
1) Mass flow rate of the cycle 
2) Number of units in parallel 
3) Shaft revolutions per minute 
4) Number of radial rotor blades 
5) Cycle conditions of section 2.2. 
2.4.4 Rotor Calculations 
The first step in determining the rotor parameters 1s to 
determine the conditions at rotor inlet (Station 3). Borrow- 
ing an analysis from the field of centrifugal compressor 
design, it can be shown that the optimum performance 1s ob- 


tained when the following relationships hold at rotor inlet. (8) 


pe= th - 2.0 


where Zi1s the number of radial rotor blades. 








ee, eee _aeee _ 
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This factor (2) is often called the slip factor. If the 
velocity at rotor exit is assumed to be entirely in the axial 
direction for reasons similar to those given for the axial 
turbine in section 2.3.4, then the slip factor may be used in 
the Euler equation to solve for the radius at rotor inlet. 
oh, = TU," 


gaa 





where U, =WRPM r,/30. 


Substituting the constants into this gives: 


+ 
r= 1510.8233\ (An)? 
3 RPM — 3 

CA 


The condition necessary to avoid negative velocities at rotor 


inlet (2) is: 


<< 
Lind 


m 


3 


biG 


ake 


> 


c. 


Combining this expression with that of the slip factor, it 1s 


possible to determine the optimum inlet flow angle, ~ 


. 
fan ot, = Voy = 93/03 
¥n3 ¥m3/ U3 


% - ArcTan (Z-2). 
(20m) 
It is now possible to determine the complete velocity triangles 
at erotor inlet, Figure 7. 


Using the above relationships and the definition for U5, 


V = ~ al 


m3 3 


Sw 
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Va = BIS OBLY (ah) s 
( 3% )° 


Ne = Vn3/ cS (ex, ) 


V93 = V3 sin (x) 


o>, 





ArcTan (V53-4, ) / ee) 


W 


3 = V43/c08(A3). 


Radial Turbine Velocity Triangles 


ROLor Pinlet 





Rotor Outlet 





Figure 7 
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Using the compressible fliow relationships, the flow area 
at inlet 1s obtained. The flow area and annulus area are 
again related by the cosine of the inlet flow angle as shown 
in Figure 2. The annulus area in this case, however, consists 
of the circumference times the blade height. From this the 
blade height can be calculated. 

At rotor outlet, the assumption is made, as stated pre- 
viously, that the tangential velocity is zero. Watanabe (9) 
states that there is a definite ratio of the mean outlet 
radius to inlet radius which will give the peak efficiency. 
This ratio is approximately 0.6. It has also been somewhat 
traditional to design for a relative outlet angle at the 
shroud of minus 60°. (2) This represents the compromise be- 
tween a low RPM machine and a low loss machine. Applying both 
of these criteria in a modified way, an outliet velocity tri- 
angle at mean, with a relative flow angle of minus 60° was 
chosen. This later proved to correlate well with previous 
work as shown in section 3.6 and avoided an iterative proce- 
dure for calculating rotor outiet conditions. 

Using these assumptions, the outlet velocity at mean, in 
terms of input variables, is 


a 
a | 2 ee PENS 
‘gown ~ 9.1347 (4h /3) (Radius Ratio) 


The value of 0.6 specified above is used as a first estimate 
for radius ratio. If the resulting outlet conditions do not 


meet the checks specified in section 2.4.9, then the ratio 1s 


-_—— 
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increased and outlet conditions recalculated. 
As 2), is minus 60". the other outlet conditions at 


Me all 3ar © 


U, =“WTRPM xr 
Ht 60 


Wan = V4,/ COs (4, ,)° 


4m 


The compressible flow relations are again used to de- 
termine the outlet flow area. In this case, since the outlet 
velocity is in only an axial direction, the flow area is 
equal to the annulus area. The outlet height may be calcula— 
ted by the same method used in section 2.3.5 to calculate 
axial blade height: 


outlet height = Area ann/2yr,_. 


Making the assumption that the outlet velocity is con- 
stant with respect to radius, and knowing the onesie hub and 
tip radii, the velocity triangles for these points are calcu- 
lated as previously shown, with the exception that the rela- 
tive outlet angle varies with the radius, or more simply 


Bp = ATAN ( a 


2.4.5 Stator Calculations 

The stator calculations are relatively simple. Watanabe 
(9) states that the clearunce between stator outlet and rotor 
inlet should be twice the passage height. If it is assuted 
that the blade or passage height is the same for both rotor 


and stutor, then the stator outlet radius is simply the radius 
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at rotor inlet plus two times the blade height. As in most 
cases this gap will not be large, and as it can be assumed 


that the moment of momentum, remains constant across the gap, 


ee aa V3 
ee Oe 


where r5/r, 1. 
This means that the velocity triungle at stator outlet is sub- 
stantially the same as at rotor inlet and will not be recalcu- 
lated. 

The radial spacing between blades or pitch, assuming that 
there are as many stator blades as rotor blades, is the circum- 
ference divided by the number of blades. The throat opening 
can also be determined in a manner similiar to that of 
section 2.3.5. Using a pitch to chord ratio of 0.6, the blade 
chord may be calculated. Using the assumption that Cos, = 
opening/spacing, the approximate radius at the stator inlet 


can be calculated as shown in Figure &. 
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Figure 8 


The inlet angle can now be determined using the law of Sines: 


Sin (6 ) >= ry pee eS eee 
“ak 

For the passage from stator inlet to outlet, it can be assumed 

that the fluid is incompressible. This is a valid assumption 

as the Mach number is low at rotor inlet. Since the mass flow 


rate must remain constant, the velocity change is proportional 


to the flow urea change. This can be expressed as 





D1 
Having determined the inlet angle,&%, the rest of the velocity 
triangle is also easily obtainable. 
2.4.6 Scroll Calculations | 
The scroll shape assumed is shown in Figure 9. It is 
assumed that the area is distributed radially ina linear 


relationship, thus: 


Area, = A (1-6_) 
S) max ir 
Where @ = radial angle in radians, 


The maximum area needed can be approximated if it is 
assumed that the tangential and meridional components of 
velocity do not vary as a function of the radial angle. The 
maximum area is then equal to the flow area at stator inlet 


which in turn is related by the cosine of the inlet angle. 


Radial Turbine Scroll Configuration 
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fail = 2 ryt cos(o ) 


where t is the blade height at rotor inlet 
2.4.7 liadial Efficiency 

As with the axial turbine, the theoretical determination 
of efficiency for a radial turbine is a much debated subject. 
Wood (4) and others have correlated efficiency with the non- 
dimensional specific speed mentioned previously, section 2.4.3; 
however, this does not give an exact enough answer for the 
purposes here. A more analytical method is needed. 

Heitt & Johnston (10) relate the performance of several 
specific turbines to several design variables; however, the 
method is rather a narrow one intended to experimentally 
evaluate performance as a function of specific variables, 
rather than provide an analytic efficiency determination. 
Rodgers (11) proposes a method of determining efficiency 
which accounts for losses caused by: 

Dise friction 

Nozzle shape 

totor friction 

kotor blade loading 
Rotor secondary flow 
liotor incidence 
Rotor clearance 


itcotor cooling 
Exhaust diffuser. 


WOAENAUAWNDN 
Ne ee Ne See Se Ne ee Nee Se” 


This method, while verified to within 2 percent, takes into 
account considerably more variables than are necessary for 
the purposes of this mathematical model. 


The method chosen wus that used by Rohlik (8) in his 








a 
study of the effects of turbine geometry on efficiency. This 
method is a compilation of several experimental correlations 
relating the efficiency to: 

1) Stator and boundary layer losses 

2) Clearance losses 

3) Rotor windage losses 

4) Rotor bouudary layer losses 

5) Exit velocity losses 

This method provided the needed precision and was based 

on the work of many acknowledged experts in the field; 
Balje (12), Heitt & Johnston (11), Wood (7), Stewart, Whitney 
& Wang (13) and Shepherd (14) to name a few. The correction 
to account for exit velocity, however, is omitted; thus the 
efficiency calculated is a total to total efficiency and is 
feripatib le with the efficiency calculated for the axial tur- 
bine. If a total to static efficiency is desired, it can be 


easily calculated from the data given for the outlet velocity 


triangles as: 


Mes = et 
5) 
Flee: aT 


285 Cea ss 


= 


where ),, = efficiency (total to total) 
N4, = efficiency (total to static) 


This conversion factor can also be applied to determine total 
to static efficiency for the axial turbine by using eee in 


lieu of pee 
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The mechanics of the method consist of relating each loss 
to an associated enthalpy change and then summing these in 
a fractional way to determine efficiency. The first loss con- 
sidered is that of rotor windage. Rotor inlet density is re- 
quired and also an absolute viscosity KO, at rotor inlet. 
The density is computed using the one-dimensional compressible 


flow relationships. 


ae 
7 eo 7-1) 


aes M2)? 
2 





Fe = P /RT, 
om 1/ (X%=1) 
OE Vein ramet (= (es1)* ee 
KT 


For Mach numbers less than one, the denominator approa- 


ches one and the density is then: 


fo=P_ (l- X-1 M 
Zoe 
RT 
O 


Sede 


Using this density, the absolute viscosity, (LO) can be 
calculated using the kinematic viscosity (YN) given as an 
input. 

A = PN. 


A Kkeynolds number for the rotor inlet is defined as: 





Substituting in the constants and the definition of rotative 


velocity: 
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Zoe 


R= 0.2094395 RPM r,°A,- 


The windage loss can now be computed as: 


, 2 
0.56 44D, 
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The clearance loss is expressed as: 


ee hip 
C O 


where 14 is defined as 


= Uy 
U= 2(.004r. ay -005r,,) j 
(Sager tees = arr aa 
e. 4 


The rotor boundary layer loss equations involve several 
assumptions (8) which go beyond the scope of this model; 
however, they are consistent with good practice and those 
made previously in this model. It is assumed that the re- 
sulting equations apply to this model. This holds true for the 
stator boundary layer loss equations also. 

The averayve of inlet and outlet blade spacing is defined 
as: 


SW ry tran tte) 


The solidity is defined as: 


- — le Or 


An ((r4)/r3)-1). 
S 
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A loss coefficient (e, ) is calculated by: 


= (0O.017¢) (1+1.9S _) 


On as OROOoZ-OnOl To) i (t5+t,) 


and the total boundary laver loss is given by: 
2 


r 4m ‘ 


e 
l- 28, J 


L = 
r 


er 
In a like manner for the rotor, a stagger angle is defined 
as: 


s Hi 
Sas (oc +5, ) /2 


where «? is given by 


x = ArcTan(Sins, 
(= > 
Ee + Co sex, j 


es 


The loss coefficient a) is calculated from 


e. = O.0076 (1+Cos <<, ) 
Cosx,-0.025 ( re 


and a total loss for the stator boundary layer calculated 
from 
2 

Ie 


a 


Ss wus 


e 
l-e. 2, py tn 
The efficiency is found from 


Nez = Abb -L /ish th +b). 


2.4.8 Miscellaneous Parameters 
Specific speed seems to be defined differently by each 


researcher. Some definitions are non-dimensional and others 








aT 
are not. Wood (7) and Kohlik (8) define their specific speed 


Similiar to that used in centrifugal pump theory as 
1 
— 2 
Ne = RPM(Q,) 
3/4 
where Q, is the volumetric flow/sec 


H is the ideal head drop across 
the turbine in feet. 


This, however, is dimensional and for that reason the 
non-dimensional specific speed recommended by Carmichael (2) 
is used instead. 


This is defined as 





The rotor outlet volumetric flow can be calculated as 
Q, = hi/P, 
where /, is obtained using the compressible 
flow relationship of the previous section. 
As a matter of interest the specific speed of Wood et al 
may be obtained from that used in the model bv using the 


following formula: 


N. (Wood) = 810.04455 N_ (Carmichael). 


The axial unit length of the rotor is a parameter which 
must be calculated in detail for the final desipn; however, in 
order to simplify the calculations in the mathematical model 


@ simpler method is used. Based on the data given in several 
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of the references, on the average, the unit axial length is 
on the order of the inlet radius less the outlet hub radius, 
thus for the purposes here, 


length = Pa-Tyane 


The machine diameters given as output are simply calcu- 
lated from the previously derived radii. 
2.4.9 Radial Limitations 

The specific speed parameter has an important influence 
on the shape of the turbine. In order that the model not 
generate dimensions which were unacceptable from a geometric 
point of view (i.e. a radical unconventional design in which 
other factors such as shroud clearance might have an unaccount-— 
ed for effect), the input parameters are first checked to de- 
termine the specific speed of the proposed turbine. If this 
specific speed is below 0.05 or above 0.25 the input para- 
meters will not produce a design which is conventionally 
feasible. 

The calculation of stress in the rotor is a complicated 
subject and one which is beyond the scope of the preliminary 
design phase. Wood (7), however, provides a means of deter- 
mining whether a design will be stress limited. Specific 
Speed is written in terms of outlet area, RPM, rotor tip speed 
and several fluid velocities. Maximum expected values for a 
stress limited case are then used to solve for a specific 


speed. This specific speed was approximately twice that 
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which would be allowed by the previous criteria and, therefore, 
the assumption is made that the radial turbines allowed by 
this model will not be stress limited. 

As in the axial turbine a continuity check must be made 
at rotor exit to assure that negative diameters are not calcu- 
lated. The rotor outlet radius is checked to insure that it 
is at least 10 percent of the inlet radius. This differs 
from the previous 50 percent used in the axial turbine; however, 
the radial rotor doesn't have to incorporate a through shaft 
and thus can be of a much smaller inner diameter. If the 
radius in question is too small, the radius ratio used in 


calculating V section 2.4.4, 1s increased by 2 percent. 


m4’ 
This has the effect of increasing the density of the fluid 
and thus decreasing the outlet annulus area required. 

The calculations described here for the radial turbine 
do not account for any shock effects caused by operation at 
a Mach number greater than 1. tven without this correlation, 
however, Rodgers (11) states that Mach number effects are 
not important until M = 1.2. Carmichael (2) recommends a 
limit of 1.1. The lower limit was used to be consistent 
with other assumptions made previously. The two velocities 
likely to be the largest, the outlet tip relative velocity, 
and the rotor inlet velocity are checked to ere that they 


are below a Mach number of 1.1. If not, the input parameters 


result in an unt'feasible turbine. 
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EVALUATION OF THE MODEL 

3.1 General 

In order to facilitate. the preliminary design of turbines 
using the mathematical model developed in the previous chapter, 
it was necessary first to program the model, and then to run 
the program over a wide range of inputs in order to verify the 
results with existing correlations; or in the case of the sim- 
pler geometric relationships, with hand calculations. the 
areas which were examined for validity and the methods used 


to establish same can be summarized as: 


Area Method 
Simple mathematical hand calculations 
relationships 
Aspect ratio selection hand calculations of 


existing turbine 


soderberg correlation comparison with 
function original function 
Axial turbine efficiency comparison with other 


published correlations 


§ 


Radial turbine efficiency comparison with other 
published correlations 


Bach of these areas will be described in more detail in the 
following sections. 

The cycle chosen for the validation was the Zener sea 
cycle as proposed by Clarence Zener CLS) eh = cycle is 
similiar to others (16) which seek to take advantage of the 
temperature differences existing in the world's oceans, and 


is representative of the type of cycle for which the mathe- 
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matical model was designed to be ead to evaluate. The de- 
tails of this cycle and of the individual parameters which 
were examined are contained. in Appendix F. It is sufficient 
to say here that over 50 different combinations of input para- 
Ineters were examined for each of two different working fluids. 
This range of input variables provided a representative 
picture of the design surface, and allowed the results to be 
broad enough to correlate with existing data over a wide 
range. 

3.2 Simple Mathematical Relationships 

This facet of the model is the predominent one, as 
explained in section 1.3, and also the easiest to verify. 

The major portion of the model is based on traditional theory 
which has been well proven in the laboratory. It is not the 
purpose here to attempt to revalidate the conclusions of re- 
searchers far more qualified than the author. The emphasis 
has been placed instead on choosing accepted relationships 
and correlating them into a usable model. 

Based on this, the validation of this segment consisted 
of hand calculating the various mathematical relationships 
and verifying the results of the computer program. To this 
end, the model has been debugged of all forseeable mechanical 
type mathematical errors. 


3.3 Aspect Katio Vorrelation 


The deeision to develop an internal choice of aspect 
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ratio based on blade bending stress for the axial turbine, 
rather than have the designer predetermine it, was based on 
primary and secondary input. factors as explained in section 
2.3.3. it was necessary then to verify that the internal 
choice was a reasonable one, consistent with existing designs. 

This was carried out by applying the correlation developed 
in section 2.3.4 to an existing design for an intermediate 
stage turbine, developed for marine propulsion use and des~ 
cribed by T. B. Hutchinson in reference 17. This design was 
for an 8 stage steam turbine ppereting at 5036 RPM and pro- 
ducing 8940 horsepower at 90 percent efficiency with an inlet 
Pe meratire and pressure of 1000°F and 225 psig respectively. 
Using the data supplied by Hutchinson in reference 17, and 
making assumptions, where necessary, of input variables which 
would optimize efficiency, the correlation of section 2.3.4 
was applied to the design and an aspect ratio was computed. 
This aspect ratio was approximately 1.4. The aspect ratio 
chosen in the actual design was approximately 1.5. 

While this simple point verification does not guarantee 
validity over the complete range of input variables, it does 
show that the relationship is reasonable. When this relation- 
ship is combined with an upper limit for aspect ratio, to take 
into account the fact that other factors than stress may limit 
the design, the resulting range of variation of aspect ratio 


(0.25 to 10.0) is consistent with existing designs (5) and is 
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assumed to provide a reasonable result. 
3.4 Soderberg Correlation 

In determining efficiency for the axial turbine, section 
2.3.5, it was necessary to formulate a function which approx- 
imated the curve of emperical data given by Soderberg (4) for 
the correlation of gas deflection angle and the primary loss 
factor. This was accomplished in a computer subprogram 
through the use of LaUrange polynomials to curve fit the 
given data to a third order approximation of the function 
(18). Without explaining the details of the theory of curve 
fitting, it is sufficient to state that within a given sub- 
interval the correlation was represented by the function 


a = x, Y, = X, Y, + X, Y, 


where Xys x, and X, are third order functions dependent upon 


be and Y.,. are values 


ie 3 


of the ordinate which bracket the given abscissa value. 


the value of the abscissa given,and Y 


The simulation of this function is shown’ in Figure 10. 
b> Axial Turbine Efficieneyv 

Although the correlation adopted for determination of 
axial total to total efficiency has been previously validated, 
it was necessary to compare the results obtained with the 
model to other correlations in order to assure. that the 
correlations provided reasonable results for the assumptions 
made in the model. This was accomplished by comparing the 


computed efficiency of several example cases for the Gener 
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cycle with the correlations presented by Craig and Cox (3). 
This comparison is shown in Figure ll. 

This comparison does not account for the differences 
between the two correlations in their accounting for secon- 
dary loss parameters; however, the results are close enough 
to show that the application of Soderberg's method is valid 
and consistent with accepted practices. 
moor Kadyval lurbine Efficiency 

It was also necessary to compare the computed radial 
turbine efficiency to other correlations in order to verify 
its application in the model. This can be more easily accom- 
plished in the radial turbine than the axial through the use 
of the non-dimensional specific speed parameter mentioned 
in section 3.4.8. Several examples, covering a wide range 
of specific speeds, were compared with the correlations of 
Wood (2) and Rohlik (8). The results of this correlation 
are shown in Figure 12. The close agreement over the wide 
range of specific speeds shows that the determination of 
radial efficiency is also valid and consistent with accepted 


practice. 
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Radial Turbine Efficiency Comparison 
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CONCLUSIONS 

The formulation of conclusions about a project of this 
nature 1s a somewhat limited undertaking. The object of the 
thesis is to compile existing data, relationships, and methods 
in the field of turbine design and to formulate and compute- 
rize a mathematical model, based on this compilation. This 
has been accomplished as shown in chapter 2. 

The validity of the mathematical model has been estab- 
lished in chapter 3 to the extent that it is valid, consistent 
with the assumptions stated in its formulation, and within 
the requirements posed by the preliminary design stage. 

The actual use of this programmed model, as explained 
in Appendix F, has shown that through the utilization of 
this program, or programs of this type, the designer may 
explore a wide range of input variables in a comparably 
short amount of time. This enables him to more fully examine 
the design spectrum and to develop a qualitatively better 


estimation of thermal cycle feasibility. 
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LECOMMENDAT IONS 

In developing the mathematical model described in chapter 
2, it became apparent that there were four areas where the 
model could be expanded which would increase its usefulness 
to the designer. the first of these is in the area of the 
flow relationships utilized. The one dimensional flow rela- 
tionships do not completely describe the flow. Expansion of 
these relationships to include shock effects would broaden 
the applicability of the model. 

The second and third areas which could be expanded are 
integrally related to one another. They entail enlarging 
the areas of the machine which are described by the model, 
perhaps including material limits and selection, disc design, 
bearing and shafting design. VUoncurrent with this would be 
the increasing number of input variables necessary to des- 
cribe the machine. This would be limited by the amount of 
detail required to satisfy a particular designer's definition 
of preliminary design. 

The last area which was not included in this model is 
that of cost estimation. Information in this area is gene- 
rally proprietory and not easily obtainable. Most of the 
information available is limited in scope and not applicable 
to the general case as described herein. This area in itself 
could be the subject of a dissertation, but certainly more 
Perk ils needed before it is acceptable for inclusion ina 


model similar to that developed in chapter 2. 
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APPENDIX A 


A.] Variable List 

The following is a list of variables as they appear in 
the main program and the subprograms. They are listed in the 
order in which they appear within the body of the programs. 
Subscripted variables refer to machine inlet and outlet, 
subscripts 1 and 2 respectively. The units of these variables 
are those listed in the Notation section, pages 7 to 9, with 
the exception of pressure, Kinematic viscosity and all ey 
which are internally converted to ine, Ne We fac and 


radians respectively. 


MAIN PiOGRAM 


Cycle and Machine Inputs 


HO] inlet stagnation enthalpy 

LO stagnation temperature 

PO stagnation pressure 

NU kinematic viscosity 

HO2 outlet stagnation enthalpy 

HO2S outlet isentropic stagnation enuthalpy 
TO2S outlet isentropic stagnation temperature 
E025 outlet isentropic stagnation pressure 

CP specific heat at constant pressure 

GAMA Pete pooLespecific heats 


CV specific heat at constant volume 





DHO 
RJ 

MM 
AXN 
RoAM 
ror MX 
ALFAMN 
UNLTX 
IePMh 
UNITR 
ZH 


WW 


ppl aight 


ttADM 


VTSMN 


ViaN 

VX 

Wasik 
TALEFNN 
BETAMN 
Bis TAMit 


ZETANL 
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universal gas constant 
stagnation enthalpy drop across machine 
universal gas constant times 778 
number of combinations of input parameters entered 
number of turbine stapes-axial 
reaction at mean 
shaft revolutions per minute-axial 
rotor inlet true flow angle at mean 
number of axial turbines in parallel on flow path 
shaft revolutions per minute-radial 
number of radiai turbines in parallel on flow path 
number of radial blades 


total mass flow rate of working fluid 


Turbine 


radius at mean 


tangential velocity component at rotor inlet, 
mean radius 


true velocity at rotor inlet, mean radius 

axlal velocity component 

relative velocity at rotor outlet, mean radius 
tangent of ALIFAMN 

relative flow angle at rotor inlet, mean radius 
relative flow angle at rotor outlet, mean radius 


primary Soderberg efficiency parumeter at rotor 
inlet 


~ 
a 
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ZETARI primary Soderberg efficiency parameter at rotor 
outlet 

CON product or radius times tangential velocity 
component 

WMN | relative velocity at rotor inlet, mean radius 

UMN tangential rotor velocity at mean radius 

FLFN function f (section 2.2.3) 

AFLO flow area 

AANN annulus area 

T blade height~axial 

KADT tip radius 

RADH hub radius 

VTTIN tangential velocity component at rotor inlet, 


tip radius 


VTHN tangential velocity component at rotor inlet, 
hub radius 


ALFATN rotor inlet true flow angle at tip radius 
ALFAHN rotor inlet true flow angle at hub radius 

VTN true velocity at rotor inlet, tip radius 

VHUN true velocity at rotor inlet, hub radius 

UTN tangential rotor velocity at tip radius 

UHN tangential rotor velocity at hub Pedi - 

WTTN relative tangential velocity component at rotor 


inlet, tip radius 


WTHN- relative tangential velocity component at rotor 
inlet, hub radius 


BETATN relative flow angle at rotor inlet, tip radius 


BE TAHN relative flow angle at rotor inlet, hub radius 





wIN 
WHN 
BETATR 
BETAHL 
WTR 
WHER 
CHD 
ASHAT 
SPN 
TBTAMN 
TBTAMK 
SPR 

TR 

TN 

DR 

DN 

RER 
REN 


LETAL2 


ZETAN2 


LETAL3 


LETAN3 


Ere 


we. 
relative velocity at rotor inlet, tip radius 
relative velocity at rotor inlet, hub radius 
relative flow angle at rotor outlet, tip radius 
relative flow angle at rotor outlet, hub radius 
relative velocity at rotor outlet, tip radius 
relative velocity at rotor outlet, hub radius 
blade chord 
aspect ratio 
stator spacing 
tangent of BETAMN 
tangent of BETAMK 
rotor spacing 
throat opening for the rotor 
throat opening for the stator 
hydraulic diameter for rotor 
hydraulic diameter for stator 
Reynolds number for rotor 
Reynolds number for stator 


soderberg efficiency parameter corrected for 
aspect ratio for rotor 


Soderberg efficiency parameter corrected for 
aspect ratio for stator 


Soderberg efficiency parameter corrected for Ak 
and Re, for rotor i 


Soderberg efficiency parameter corrected for Ak 
and ke, for stator 


staye efficiency 





LX 
NZX 
EPPX 


AMN 


AMR 


hEAH2 
ESTMX 
DIAM1 
DLAM2 
UNLEN 


G 
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number blades-decimal number 
number blades-whole number 
machine efficiency 


true Mach number at rotor inlet, hub radius, last 
stage 


relative Mach number at rotor outlet, tip radius, 
last stage 


reaction at last stage, hub radius 
centrifugal stress at blade root, last stage 
inlet diameter 

outlet diameter 

machine length 


radians to degrees conversion 


itadial Turbine 


Z 
HAD3 
ALIA3 
v3 
V3 
AM3 
FUN3 
AFLO3 
TR 
hAD2 


KirAT 


slip tactor 

radius at rotor inlet 

true inlet flow angle, rotor inlet 
meridional velocity component at rotor inlet 
velocity at rotor inlet 

Mach number at rotor inlet 

function f ‘at rotor inlet (section 2.2.3) 
flow area at rotor inlet 

stator blade height 

radius at stator exit 


rudius ratio (section 2.4.4) 





V4M 
V4H 
HAD4M 
XM4M 
FUN4 
AANN4 
T4 
WAD4H 
RAD4T 
BB 
DEN4 
Q 
SPSP 
SPI 
UP 
CK 
RADI 
SALFAI 
CALFAI 
ASC 
V1 
VMI 
VT) 
ALFA] 


wt3 


17 
velocity at rotor outlet, mean radius 
velocity at rotor outlet, hub radius 
mean radius at rotor outlet 
Mach number at rotor outlet 
function f at rotor outlet (section 2.2.3) 
annulus area at rotor outlet 
blade height at rotor outlet 
hub radius at rotor outlet 
tip radius at rotor outlet 
exponential constant 
fluid static density at rotor outlet 
volumetric flow rate at rotor outlet 
non-dimensional specific speed 
stator outlet spacing 
stator throat opening 
stator chord 
stator inlet radius 
sine of ALFA] 
cosine of ALFA] 
maximum scroll] area 
velocity at stator inlet 
meridional velocity component at stator inlet 
tangential 
true flow angle at stator inlet 


tangential velocity component at rotor inlet 





U3 
BETA3 
W3 

U4a1 
U4H 
U4T 

bi TA4M 
BETA4H 
BETA4T 
W4h 
W4H 
W4T 
XM4 
KhE3 
DEN3 
XLM 
ALC 
SPM 
SLGMAR 
ket 

XLR 
ALFAO 
ALFAST 
ES 


XLS 
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tangential rotor velocity at inlet 
relative flow angle at rotor inlet 
relative velocity at rotor inlet 
tangential rotor velocity at outlet, mean radius 
tangential rotor velocity at outlet, hub radius 
tangential rotor velocity at outlet, tip radius 
relative flow angle at rotor outlet, mean radius 
relative flow angle at rotor outlet, hub radius 
relative flow angle at rotor outlet, tip radius 
relative velocity at rotor outlet, mean radius 
relative velocity at rotor outlet, hub radius 
relative velocity at rotor outlet, tip radius 
relative Mach number at rotor outlet, tip radius 
Reynolds number at rotor inlet 
fluid static density at rotor inlet 
windage loss 
clearance loss 
mean rotor spacing 
solidity 
rotor boundary layer loss coefficient 
rotor boundary layer loss 
stator angle (section De lene) 
stator stagger angle 
stator boundary layer loss coefficient 


stator boundary layer loss 





HEE 
UNLENK 
DIA3 
DIA4T 
DIA4M 
DIAl 
DIA4N 


DIA2 
SUBPRKOGHAM 


Xl 
X 
XX 
YY 
ALL 
XL2 


XL3 


SUBPicOGRAM 


GAMA 
ht 

CP 
TO 


V 


(oe, 
machine efficiency 
machine length 
diameter at rotor inlet 
tip dianeter at rotor outlet 
hub diameter at rotor outlet 
diameter at stator inlet 
mean diameter at rotor outlet 


diameter at stator outlet 


SODCOR 


gas deflection angle, radians 

gas deflection angle, degrees 
abscissa for Soderberg correlation 
ordinate 

first LaGrange Polynomial 

second LaGrange Polynomial 

third LaGrange Polynomial 


Soderberg gas deflection loss parameter 


PLOREL 


ratio of specific heats 

universal gas constant 

specific heat at constant pressure 
Stagnation temperature 


velocity 





XM 
FLEIEN 
4 
EXPR 
CON 


FUN 


Mach number 

function f (section 2.2.3) 
static temperature 
internal constant 

internal constant 


internal constant 


80 





81 


APPENDIX B 


B.1  =Main Prograin 


A complete listing of the main program is contuined on 


the following pages. Appendix E contains detailed instruc- 


tions for input-output and running the program. The inathe- 


Inatical model for the program is developed in chapter 2. 


The logic sequence used in programming the model is outlined 


below: 


Oecenanuwbwhd He 
Nee” Snel” “iene” “ieee” “emt” “inet” “ieee” “ene” “eet” “tenet” 


~~ 


a 
b= 
— 


12) 
13) 


14) 
15 ) 
16) 
17) 
18) 


19) 
PO) 
21 


Definition of real variables and common variables 
Dimensioning of subscripted variables 

Input format specification 

Gutput format specification 

Readin input data 

leadout input data for checking 

Convert input data to standard units 

Calculate fluid properties 

kkeadout fluid properties for checking 

Commence do-loop for calculations on each set of 
turbine input parameters 

Calculate axial turbine mean parameters 

Calculate axial first and last stage parameters on 
internal do-loop 

If design is rejected on specified checks, axial cal- 
culations are dropped and radial are commenced. 
Calculate axial machine parameters 

Convert axiul variables to output units 

leadout output data for axial inachine 

Calculate radial parameters 

lf design is rejected on specified checks, radial 
calculations are dropped and next set of axial 
calculations commenced 

Convert radial variables to output units 

kteadout output data for radial machine 

Pommence ca leuvlaeilonsmrorgnex. SeLmofeaxtaly input 
parameters 
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APPENDIX C 
CC Subprogram SODCOR 
This program computes ‘an efficiency parameter from 
given data points using LaGrange Polynomials in a curve- 
fitting logic. <A complete listing of the subprogram is 
contained on the following page. The mathematical model 
is explained in section 3.4. The logic sequence used in 


programming the model is outlined below: 


1) Transfer of input variables 

2) Definition of curve data points 

SCONCE StoOnTor input data to ™standard units 

4) Classification of input data as to sub-division 
of curve for curve-fitting purposes 

5) Definition of sub-division curve-fitting data 

6) Calculation of LaGrange Polynomials 

7) Calculation of output variable 

8) Transfer of output variable 
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APPENDIX D 

D Subprogram FLOREL 

This subprogram computes a flow function and Mach number 
using the one-dimensional flow functions. The mathematical 
model is explained in section 2.2.3. A complete listing of 
the program is contained on the following page. The logic 
sequence used in programming the model is outlined below. 

1) Transfer of input variables 

2) ‘Transfer of common variables 

3) Calculation of Mach number 


4) Calculation of flow function 
5) Transfer of output variables 
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APPENDIX E 

k.l Specific Instructions for Use 

This section is intended to provide detailed instruc-— 
tions for the use of the computer programs described in the 
previous appendices including input-output interpretation. 
The specific program control cards described are required 
on the Interdata 7O computer. Check the user's munual be-~ 
fore running on other machines as these control cards may 
differ. 
k.e2 Loading Sequence 

The sequence for loading the main and subprograms with 


their required control cards is as follows: 


1) Job card 

2) Language card 

3) Listing designation card 

4) Main program 

5) Subprogram SODCOR 

6) Subprogram FLOREL 

7) Execution card 

8) Common length designation card 
9) Data cards 
10) Termination card 


B.3 Job Vontrol Cards 

- The JOB card, faneuatc card, listing designation card 
and JOB’ termination cards are all of the standard format, 
eet any job run on the Interdata 70 system. The ex- 
ecution card must carry a numeral one in column seventeen in 
addition to the standard format. ‘The common length desipna-— 
tion card specifies the space in blank common to be used in 


the job. The format for this card is *BC OQOOC commencing 
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in column one. 
£.4 Data Cards 
BE.4.1 General 

The Data cards necessary for use of the program can be 
grouped into three categories. 

1) Cycle Data 

2) Number of Design Sets 

3) Turbine Design input sets 
They must be loaded in the above order. Card groups 1 and 2 
are only required once per job, while card group 3 is re- 
quired for each individual design set examined. 
BE.4.2 Cycle Data 

Three cards are required to specify the data for a 
given thermal cycle. The context of each is given below. 
No information can be omitted, and they must be loaded in 


& 


the order described. 


SPACKS INFORMATION UNITS FORMAT 
CARD ONE 

1-10 inlet total enthalpy BLU/1bm F10.5 
20-30 inlet total temperature OR H1LO% 
40-50 inlet total pressure psia EVO=> 
60-70 inlet kinematic viscosity et~/hr PO 5 
CARD TWO 

1~10 outlet total enthalpy BTU/1bm PLO o 
20-30 outlet total temperature OR Cs S 
40-50 outlet total pressure psia F10.5 
60-70 © outlet kinematic viscosity ey re F10.5 
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SPACES INFORMAT LON UNITS FORMAT 


CARD THREE 


1-10 outlet isentropic total BTU/1bm F10.5 
enthalpy 

20-30 outlet isentropic total me F10.5 
temperature 

40-50 outlet isentropic total psia F10.5 


pressure 
&B.4.3 Number of Design Sets 

One card 1s necessary to specify the number of sets of 
turbine parameters entered. A set consists of one combina-— 
tion of axial turbine input parameters, and one combination 
of radial turbine input parameters, with both combinations 
sharing the same total mass flow rate (section 2.2.2). The 
program will accommodate up to 99 separate sets of input 
data. The number of sets must be specified in [2 format on 
spaces one and two of the card with the other spaces left 
Tanke This card may not be omitted. 
E.4.4 Turbine Design Input Sets 

A turbine design input set consists of a combination of 
axial input parameters, section 2.3.3, and a combination of 
radial input parameters, section 2.4.3, which share the same 
total@mass filliow Tate, section 2.2.24. Three cards! are 
necessary for each set and no information or card may be 
omitted. they must be loaded in the order in which they are 
described. In the event that a run is not to be designated 


by name, a blank card must be included in place of card one. 





SPACES 
CARD ONE 


6-70 


INFORMATION UNITS 


Any desired run designa- 
tion intormation 


CARD TWO-Axial Turbine 


11-15 
21-25 
31-40 
51-55 


61-65 


Number of stages 
Inlet reaction at midblade 
Shaft RPM Rev/min 


Rotor inlet flow angle Degrees 
at midblade 


Number of machines op- 
erating in parallel on 
total flow path (may be one) 


CARD THREE-Radial Inflow Turbine 


11-20 
31-35 


41-45 
51-60 


Bes Cutput 


Shaft RPM Kkev/min 


Number of nachines op- 
erating in parallel on 
total flow path (may be one) 


Number of rotor blades 


Total mass flow rate for lbm/sec 
all machines—-i1.e. sum of 

all parallel flow paths. 

see section 2.2.2. ee 

must be the same for Axial 

and Radial Design Set. 


Po. beecreneral 


The cycle data will be printed out once for each 


in the same format 1t was entered in with, 


accuracy. 
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FORMAT 


LOSES) 
1805 ore 
me OI GIS. 
be. 2 


F5.0 


F10.0 
F5.0 


F'5 .0 
res LON 


job, 
as a check for 


The universal gas constant, in units of ft lbf/ 


lbm°k, the specific heat at constant pressure, and the ratio 


of specific heats will also be printed out as a check. 


ror 
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each turbine design data set, the run designation information, 
the assumed input variables, and the preliminary jestne data 
will be printed out. If the assumed input variables result 
in a design which doesn't pass the axial or radial feasibi- 
lity checks, sections 2.3.8 and 2.4.9, the design data will 
be omitted and an unfeasibility message printed, section 


— 


hoe Be 
BE.5.2 Preliminary Design Data 

The format of the preliminary design data is self- 
explanatory. The output is in units of feet for all dimen- 
sions, degrees for all angles, and feet per second for all 
velocities. stress 1S given in psi and tpi. The abbrevia- 
tions used in describing the velocity triangles are as 
follows: 
Axial Turbine 


UHN, UMN, UTN, - rotor tangential velocity at hub, mean and 
tip radii 


VTHN,VITMN,VTUN,- fluid tangential velocity component at 
rotor inlet at hub, mean and tip radii 


VX ~ fluid axial velocity component 


VHN, VMN, VTN, - fluid velocity at rotor inlet at hub, mean 
and tip radii 


WHN, WMN, WIN, —- fluid relative velocity at rotor inlet at 
~ hub, mean and tip radii 


WHR, MMR, WTR, -—- fluid relative velocity at rotor outlet, hub, 
mean and tip radil 


ALFAHN, ALFAMN,ALFATN - rotor inlet flow angle at hub, mean and 
tip radii 
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BETAHN, BETAMN,BETATN, - relative rotor inlet flow angle at 
hub, mean and tip radii 


BETAHK,BETAMR,BETATR, - relative rotor outlet flow angle at 
hub, mean and tip radii 


Radial Turbine 


V1, V3 — fluid velocity at station 1 or 3 

celine V 13, - fluid tangential velocity component at 
station 1 or 3 

vile VM, - fluid meridional velocity component at 
station 1 or 3 

U3 - rotor tangential velocity at inlet 

W3 ~ fluid relative velocity at rotor inlet 

ALFA1,ALFA3, - inlet flow angle at station 1 or 3 

BETA3 - relative rotor inlet flow angle 

U4H,U4i,U4T, - rotor outlet tangential velocity at hub, 


mean and tip radii 


V4H,V4N,V4T, - rotor outlet fluid velocity at hub, mean 
and tip radil 


W4H,W4N,W4T, - rotor outlet fluid relative velocity at 
: hub, mean and tip radii 


BETA4H, BETA4M,BETA4T, - rotor outlet fluid relative flow 
angle at hub, mean and tip radii 


B.5.3  Unfeasibility Conditions 

If a given set of turbine design inputs does not meet 
Mite fieash bility checlds delscribed in sections 2.3.8 and 2.4.9, 
a design unfeasibile message wili be printed. As a reference, 
the checks associated with each messave are summarized below. 


MESSAGE INTERPRETATION 


Axial Jurbine 


Gas Deflection too big Bither %, or SP +5) is greater 
than 140° 
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MESSAGE INTERPRETATION 

Hub Kadius too small ry < 0-5 ry 

Hub iteal Mach Number Last stage hub true Mach 

too big number greater than 0.99 

Tip Real Mach-Number Last stage tip relative Mach 

too big number greater than 0.99 

Hub reaction too small Last stage hub reaction less 
than 0.01 


Kadial Turbine 


Specific Speed to Radical Machine specific speed less 
than 0.05 or greater than 
O25. 

Inlet Mach number too Mach number at rotor inlet 

high greater than 1.1. 

Outlet Relative Mach Itelative Mach number at rotor 

number too high tip greater than 1.1 


k.5.4 Compilation Errors 

It is possible, given certain sets of axial turbine in- 
put data, to cause a computer run time error (negative square 
root and overflow). This condition sometimes occurs during 
the calculation of aspect ratio. It has been forseen in pro- 
gramming the model and will be internally corrected if it 
occurs. The data output is still valid if this error occurs 
even though an error message is printed out. Other run time 
errors or errors in computing the radial turbine data should 


be cause for re-examination of the input data provided. 
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APPENDIX F 
EVALUATION OF THE ZENER CYCLE 

F.l General 

The Zener Cycle, as described by Clarence Zener (15), 
1s based upon the thermocline existing in the oceans of the 
world. According to Zener, it should be possible to take 
advantage of this temperature difference, (20°C is used as 
an example) to power a thermal cycle which could in turn 
produce useful energy. This supposition is based upon the 
assumption that technology could produce the necessary 
equipment. Assuming that this is possible, Zener estimates 
that "the tropical oceans in the year 2000 could supply the 
whole world with energy at a per capita rate of consumption 
equal to the U.S. per capita rate in 1970, and suffer only a 
one-degree C drop in temperature". This is estimated to be 
about 60 biliion kilowatts. Zener proposes a thermal cycle 
using emon ie as the working fluid and operating between ee 
and 25°C. His estimated plant size is about 8000 cubic feet. 

The cycle consists of taking high pressure liquid at 
BS Cs expanding 1t through a turbine, and condensing the 
vapor at 5°c. The fluid would then be pressurized and 
heated to 25°C to begin the cycle again. Sea water would 
be used to both heat and cool the fluid through the use of 


heat exchangers. Ammonia is suggested as a working medium. 
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F.2 Working Fluid 

Without exploring the heat exchanger problems involved 
in detail, the feasibility of the turbines needed for such 
a cycle was evaluated. The calculations were based upon a 
plant output of approximately 100,000 kilowatts with an 
assumed turbine efficiency of 90 percent. Based upon this 
the flow rate for a specific fluid was calculated. Ammonia 
and Freon 2] were chosen as the fluids; ammonia because it 
was proposed by Zener, and Freon 21] because it is typical of 
"the recently developed refrigerating fluids" mentioned by 


Zener. The cycle operating points are summarized as follows, 


along with the necessary mass flow rates. 


Point Enthal py Tem Pressure 
(BLU/1bm) (On) (psia) 

AMMONIA 

] 628.74 528 | 124.4 

2 612.575 510 89.19 

PAN 610.779 510 599 LS) 
FREON 2] 

] 27 306 523 Pee 

2 24 25 510 5 Oo 

25 ts. IGS 510 5 230 


MASS FLOW RATE 
Ammonia — 6050 !tbm/sec 


Freon 2) — 30220 lbm/sec 
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F.3  kange of Input 

The inputs chosen for the radial turbine were based 
upon a back calculation from specific speed to RPM and number 
of units in parallel. The RPM was quantified as being 3600/N 
where N is an integer number. This was done because it is 
assumed that the machine would be used to produce 60 hertz 
power without benefit of a reduction gear. The number of 
units was then specified to provide a variety of machines 
covering the range of specific speed recognized as feasible, 
section 2.4.9. Twelve rotor blades were specified for the 
ammonia trials and 15 for Freon 21. This brackets the usual 
design range. 

The axial turbine inputs were based on combinations of 
parameters used to achieve high efficiency. A large slow 
machine was favored with high reactions and few stages. 

This was not immediately apparent at the outset, but after 
several runs with no combinations yielding a feasible turbine, 
it was indicated. 

A summary of design inputs for one run is shown below. 
The feasibility column indicates whether the design meets 
the feasibility criteria of the model, sections 2.3.8 and 
24.9. 

Axial Turbine Inputs-Ammonia 


# Stages KKeaction RPM a Units Feasible 


2a a, 900 60 4 Yes 
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# Stages keaction RPM an Units Feasible 
10 » - “aco 60 4 No 
5 aD 720 60 4 Yes 
2 a, 720 60 4 Yes 
10 ms, 900 60 2 No 
y) a, 900 60 2 No 
2 “eS, 900 60 2 Yes 
10 a, 720 60 2 No 
) 45 720 60 2: No 
2 a6, &2O _ 60 2 Yes 
5 a, 900 60 4 Yes 
10 25 900 60 4 No 
2 a5) i2o 60 10 Yes 
5 Pap, 720 60 10 Yes 
10 so 720 60 10 Yes 
2 oo 900 60 10 Yes 
») 55) 900 60 10 Yes 
10 , 900 60 10 No 
5 <5 900 60 4 Yes 
) Ae 900 60 4 Yes 
> of 900 60 4 Yes 
y) aS 900 60 4 Yes 
y) al 900 70 4 No 
5 6 900 70 4 No 
y) il 900 70 4 Yes 
) 0 900 70 4 Yes 
5 aS 900 80 4 No 
5 .6 900 80 4 No 
5 ae | 900 80 4 No 
5 545) 900 80 4 No 
20 4 720 60 20 No 
20 ee 600 60 20 No 
10 we 600 60 20 Yes 
20 2 600 50 20 No 
10 ns 600 60 20 Yes 
20 = 600 60 20 No 
20 nS 600 50 20 Yes 
20 . 600 40 20 Yes 
20 a 600 30 2 Yes 
2 OF a3 600 20 20 Yes 
Axial Turbine Inputs-Freon 2] 
2 aw 600 60 20 Yes 
5 at | 600 60 20 No 
10 ot 600 60 20 No 
2 eft 20 60 20 No 
y) suf (e410) 60 20 No 
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it Stages keaction RPM ae Units Feasible 
10 a" 20 60 a No 
2 of 7, 24@) 70 20 No 
9 a7 720 rae 20 No 
10 Ags 720 AO 20 No 
2 ae 720 8O 20 No 
> 5 if 720 80 20 No 
10 ai 720 80 Zo No 
2 a) 720 60 20 No 
2 ye 720 . 60 20 No 
2 -4 720 60 20 No 
2 iG 720 60 20 No 
2 a 720 60 20 No 
2 co 720 60 20 Yes 
2 oY fF 20 60 20 No 
2 ai 600 70 ZO No 
») ni, 600 70 20 No 
10 A 600 7O XG, No 
2 Ave 600 SO 20 No 
9 aif 600 80 20 No 
10 ok 600 80 20 No 

Radial Turbine Inputs—Ammonia 

RPM Units # Blades Feasible 
3600 13 P22 Yes 
3600 2 12 Yes 
3600 1l ez Yes 
3600 10 12 Yes 
3600 9 12 ) Yes 
3600 8 12 Yes 
3600 i 12 Yes 
3600 6 12 Yes 
3600 5 12 Yes 
3600 4 2 No 
3600 14 12 Yes 
3600 | 15 12 Yes 
3600 16 We Yes 
3600 © Le 12 Yes 
3600 18 eZ Yes 
3600 19 12 Yes 
3600 20 12 Yes 
3600 AMT 12 Yes 
3600 oe 12 Yes 
3600 fa 2 Yes 
3600 24 ee Yes 


1800 lf l2 No 
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f.4 Sample Output 


which could be called feasible. 


Units 


FMR AWA WwW 


Radial Turbine Inputs-Freon Pat 


250 

300 
20 
1000 
1250 
1500 
1750 
2000 
7a Paps © 
2500 
ANAS, 
3000 
3250 
3500 
oi 28 
4000 
4250 
4500 
4750 
5000 
5250 
3300 
2150 
6000 
6250 


# Blades 


Me 


15 
15 
15 
15 
T2 
15 
15 
15 
15 
Ws 
15 
15 
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Feasible 


Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 


No 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 
Yes 


The output for both working fluids produced many designs 


A representative sample of 


axial and radial designs for each fluid is included on the 


following pages. 


These designs are not necessarily optimums, 


as the object of the runs was to test the program rather 
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than optimize a design. Data from these runs was used to 
substantiate the evaluation of the model, sections 3.5 and 
3.6. These examples do, however, illustrate several points 
about the Zener Cycle and will be the basis for the next 
section. 

Axial Turbine-Ammonia 


Assumed 5. Staves 0.50 Reaction at mean 720. RPM 
Assumed 60.00 ALFAMN 4. Units 6050.00 i 


Last stage tip rel Mach no. 0.306 Last stage hub real Mach 
no. O.314 
First stage tip diam. 8.51 Blade height 0.9599 # blades 50. 
| tao. Sooo 
Last stage tip diam. 8.84 Blade height 1.2949 # blades 38. 
Bit. 6.9593 
Machine efficiency 0.9594 
Max cent. stress 3003.50 
Unit length 2.492 
Aspect ratio, first stage 2.17 Last stage 2.34 


Velocity Triangles 

First Stage 

Hub 

UMN a45e5m VIHN 325.9 VX 164.2 VEN 365.0 WHNS@181.7 
Whit 29757" ALMAHN 63.3 BETNHN 25.3 BETAUHR =56.5 
Mean 

UMN 284. NW CUSN Ge? ney ee LOee ee VitN@ > 2 o> WIN el o4c 2 
Wh 328.5 ALFAMN 60.0 BETAMN 0.0 BETAMt -60.0 

Tip 

o Neee ome ¥ LEN 2528 a VS Lot Ze VIN@ OO Siew iN 177-9 
WTit 360.3 ALFATN 56.9 BETATN -22.6 BETATR -62.9 
Last Stage 

Hub 

UIUIN 235.6 VTUN 343.4 VX 164.2 VHN 380.6 WHN 196.4 
Witieue one camer AlINRI64 245 BETAHN 9323 Soi TAH =5'5 .1 
Mean 

OMN, 264e 0m LINING S4.5 VX 164.2) YMNEO282> WMN 164.2 


oa! 


WMt 328.5 AL#rAMNN 60.0 BLTAMN 0.0 BETAMIC -60.0 

Maas 

Dio temeev LIN 242.5 VX 164.255VIN 293 18 WIN 187.5 
5 ALFATN 55.9 BETATN -28.9 BETATR -63.8 


Wilita 7 le 
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Axial Turbine—Freon 21 


Assumed 2. Stages 0.70 iteaction at mean 720. hPM 
Assumed 60.00 ALFAMN 20. Units 30220.00 m 


Last stage tip real Mach no. 0.741 Last stage hub real 
Mach no. 0.542 

First stage tip diam. 8.87 Blade height 1.9793 # blades 84. 
tel sO ¢ lie 

Last stage tip diam. 9.67 Blade height 2.7744 # blades 23. 
pmo 0) Sire I 

Machine efficiency 0.9706 

Max cent. stress 5878.27 

Unit length 0.802 

Aspect ratio, first stage 10.00 Last stage 4.59 


Velocity Triangles 

First Stage 

Hub 

Shieh ee  VINNe2ZIS.7 = VX 90.0 VHN 23695" WHN’ 96.0 

Wilk 205.9 ALFAHNN 67.6 BETANHN 20.4 BETAHR -64.1 

Mean 

DiNmeco9soe VINN 155.9" VX 90.0 VMN 180.0 WMN 137.5 

WMR 275.0 ALFAMN 60.0 BETAMN -49.1 BETAMR -70.9 

Tip 

wee eee LIN TOP VA 9OTO VIN 1502909 WIN 231.5 

WTR 346.3 ALFATN 53.4 BETATN -67.1 BLTATR -74.9 

Last Stage 

Hub 

Dail Sowee LEN 2600,9euev% 90.0 —VIIN 276.0 WHN 338.8 

WHR 179.4 ALFAHN 71.0 BETAHN 49.6 BETAHR -59.9 

Mean . 

DMN. Ome IMN 155.95 VR 90.08 VMN 18050" WN 137.5 
O 


WMR 275.0 ALFAMN 60.0 BETAMN -49.1 BETAMK -70.9 
Tip 
UTIN 364.4 VITN 111.2 VX 90.0 VTN 143.0 WIN 268.8 


4 


peice 2 ALFATN 51.0 BLTATN -70.4 BETATR -76.1 





Radial tinflow Turbine-Ammonia 
Assumed 3600 KPm 17. Units 12. Blades 6050.00 WW 


Nozzle inlet diam. 6.99 Blade height 0.2338 
Nozzle outlet diam. 4.63 

Rotor inlet diam. 3.70 Outlet tip diam. 2.91 
Rotor hub diam. 1.52 

Unit length 1.09 

Max scroll area 3.972 

Specific speed 0.127727 

Efficiency 0.9398 


Velocity Triangles 

Nozzle Inlet 

Wao ol.G VI) 229.8 VMI 260.4 ALFA] 39.3 

Rotor Inlet 

W3)6096.8 VI3 580.7 ViB 3604.9 V3 685.38 W3 382.9 
ALFA3 57.9 BETA3 -17.7 

motor Outlet 


Hub 

WA oe aoe V4n 24174) Wal 375.3 BETA4H =50.0 
Mean 

U4hM 418.1 V4M 241.4 w4M 482.8 BETA4M ~—60.0 
Tip 


U4T 548.8 V4T 241.4 W4T 599.5 BETA4T -66.3 





Kkadial Inflow Turbine-freon 21 
Assumed 3600 RPM 1250. Units 


Nozzle inlet diam. 2.75 Blade 
Nozzle OULletedian.. 2.04 

howe anletvewdiram. 1.62 Uutlet 
hRoOLOm ilies di am. O- G2 

Unit length 0.45 

Max scroll area 0.62% 

Specific speed 0.113886 
Efficiency 0.9345 


Velocity Triangles 

Nozzle Inlet 

Vay 3.4 VIL 124.8 “Y¥M1 120.4 
Rotor Inlet 

Uae O> .Smev 13s 26570 VM3 126.1 
ARMAS O42 BETAS —~17.7 

iO GO e011 tae. 

Hub 


15. Blades 30220.00 WW 


height 0.1046 


ta pedilame 


eee 3 


ALFA] 46.0 


V3 294.3 


W3 134.4 


U4li 135.1 V4H 105.9 W4H 171.7 BETA4H ~51.9 


Mean 


U4M 183.5 V4mM 105.9 W4mM 211.8 BETA4M -60.0 


Tp 


U4T 231.8 V4T 105.9 W4T 254.9 BETA4T ~65.4 


He 
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F.5 Comment _on the Zener Cycle 
Using the sample output shown in the previous sections, 
it 1s possible to calculate a volume required for each machine. 
Multiplying this by the number of units gives a volume to 
produce 100,000 kilowatts of power. This approximate volume 


3 3 


is 600 ft” for Axial with Ammonia, 1200 ft~ for Axial with 


3 for Radial with Ammonia, 3700 ft? f@. T° 


Freon 21, 700 ft 
Radial with Freon 21. This is interesting in that it shows 
ammonia to be the more advantageous of the two fluids in 
terms of size. This would be expected, however, in that for 
the cycle points defined, ammonia has a greater enthalpy 
change per pound than Freon. Freon has the advantage of 
being safer to work with, however. 

A second fact 1s apparent in the volume consideration 
in that the estimated volume was in the order of one to two 
thousand cubic feet. Zener estimated a turbine module to be 
8 x 8 x 40 feet or about 2560 cubic feet. This shows that 
Zener's estimation was reasonable for plants of 100,000 
Kilowatt output. Assuming this size plant, then to supply 
the output of 60 billion kilowatts mentioned in section F.1, 
600,000 individual plants would be required. This also 
assumes that the heat exchangers are of comparable size. 

Cowparisons could be made for the amount of working 
medium required, amount of cooling water and others. That 


is not really the point of the evaluation, however. Ixami- 
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nation of the sample designs does show, for at least the 
turbines, that the proposal is at least within the realm of 


possibility if not reasonability. 
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